Scattering functions of knotted ring polymers 
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We discuss the scattering function of a Gaussian random polygon with N nodes under a given 
topological constraint through simulation. We obtain the Kratky plot of a Gaussian polygon of 
N — 200 having a fixed knot for some different knots such as the trivial, trefoil and figure-eight 
knots. We find that some characteristic properties of the different Kratky plots are consistent with 
the distinct values of the mean square radius of gyration for Gaussian polygons with the different 
knots. 
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I. INTRODUCTION 

Ring polymers have attracted much interest in polymer 
physics, and various properties have been studied both 
theoretically and experimentallyi 2 * 3 ^*^ For the Gaus- 
sian random polygon the analytic expression of the static 
structure factor was obtained by Casassai The Kratky 
plot of the scattering function is compared with that of 
star polymers with four or five armsi^ The scattering 
data of cyclic polystyrene in deuteriated cyclohexane was 
obtained by the SANS experiment^ 

Recently statistical properties of ring poly- 
mers under topological constraints have been in- 
vestigated extensively mainly through computer 
simulationi 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 It is first conjectured by 
des Cloizeaux that a topological constraint should 
lead to an effective repulsion among segments of ring 
polymers^. The conjecture is supported by the numerical 
observations that the mean square radius of gyration 
of ring polymers with a fixed knot is larger than that 
of no topological constraint i 8 : 10 : 11 ' 12 : 13 : 14 : 15 In fact, the 
topological swelling is observed particularly for ring 
polymers with small or zero excluded volume^ 1 ^ 

In this paper we discuss the scattering function of ring 
polymers under a topological constraint. It should be 
fundamental for studying ring polymers in scattering ex- 
periments. We consider ring polymers in solution at the 
9 temperature, and they are modeled by random poly- 
gons. Here, random polygons have no excluded volume, 
i.e. they have no thickness. We shall evaluate the radial 
distribution function in simulation, and then take the 
Fourier transformation. We shall show the Kratky plot 
of the scattering function of a random polygon having 
some fixed knot type. 

II. SIMULATION METHODS 

Making use of the conditional probability 
distribution^, we have systematically constructed 



10 5 samples of the Gaussian random polygon with 200 
nodes. We have calculated two knot invariants, A.k(— 1) 
and V2(K), to each of the 10 5 configurations, and effec- 
tively classified them into different topological classes. 
Here the symbol Ajf(-l) denotes the determinant of 
a knot K, which is given by the Alexander polynomial 
A(i) evaluated at t = —1. The symbol 1>2(K) is the 
Vassiliev invariant of the second degree. We select 
such polygons that have the same set of values of the 
two knot invariants. Here we calculate V2(K) by the 
algorithmic. The two knot invariants are practically 
useful for computer simulation of random polygons with 
a large number of polygonal nodesi&. 

We consider four different topological classes: the triv- 
ial knot (0), the trefoil knot (3i), the figure-eight knot 
(4i) and the other knots (other), in the paper. We de- 
note by "all" such polygons that have no topological con- 
straint. 

Let us denote by (R G k ) the mean square radius of 
gyration for such Gaussian polygons that have N nodes 
and a given topological constraint K . The estimates of 
(Rq k ) for N — 200 are given in Tablc[T|for several knots. 
For different numbers of N, they have been obtained in 
Ref. UJJ- The fraction of Gaussian polygons with a given 
knot has also been evaluated for some knots^ 



K (Rg,k) errors 

18.033 0.082 

3i 16.208 0.120 

4i 15.043 0.250 

other 13.459 0.118 

all 16.674 0.060 



TABLE I: Mean square radius of gyration (Rq k ) of a Gaus- 
sian random polygon with a topological condition K. Here 
N = 200. 
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All Knot 
Trivial Knot 
Trefoil Knot 
Figure Eight Knot 
the Other Knot 




FIG. 1: The probability distribution Anr 2 gK(r)Ar/N versus 
the distance r. Here N = 200 and Ar = 0.1. The plots of 
the five topological conditions, all, 0, 3i, 4i and other, are 
represented by crosses, tilted crosses, double crosses, open 
squares and open circles, respectively. 



III. RADIAL DISTRIBUTION FUNCTIONS 



K 7-pcak Rg,K 

4.25 4.246 

3i 4.05 4.026 

4i 3.85 3.879 

other 3.65 3.669 

all 4.05 4.083 

TABLE II: Peak position rp ea it of the probability distribution 
Anr 2 gx(r) Ar/N for a topological condition K. The estimates 
of r P eak may have errors of order 0.1 at most. 
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Let us define the segment pair correlation function of 
a random polygon with N nodes by22> 



1 N 

m.n— 1 



(i) 



Here R m denotes the position vector of the mth node 
for to = 1,2,...,N. Due to spherical symmetry, the 
function g(r) depends only on the distance, r = \r\, and 
we denote it by g(r). We call it the radial distribution 
function. Here we note that Airr 2 g{r)Ar/N gives the 
probability of other segments appearing in a spherical 
shell from radius r to r + Ar centered at a given segment. 

For a Gaussian polygon under a topological condi- 
tion K, we denote by gxif) and guir) the pair corre- 
lation function and the radial distribution function, re- 
spectively. The graphs of the probability distribution 
4/T r 2 gK{r) Ar/N are plotted against r in Fig. ^ They 
are consistent with a preliminary resulliSi. 

We now generalize the recent observation that for ran- 
dom polygons the peak position of the probability dis- 
tribution 4:TTr 2 g(r)A/N should be given by the gyration 
radius. 2 » Let us denote by Rq,k the square root of the 

mean-square radius of gyration, i.e. Rq.k = J {R% k) • 

The estimates of the peak position of the probability dis- 
tribution ^Trr 2 gx{r)Ar /N for the five topological con- 
ditions are listed in Table |H] together with those of 
the gyration radius, Rg,k- They are almost identical 
up to numerical errors. Thus, we have r pca k = Rg,k 
within errors. This is characteristic to ring polymers. In 
fact, for a Gaussian linear chain, the peak position of 
the probability distribution Airr 2 g(r)Ar /N is located at 
r w 0.7ARa i in , where Rg,Uu denotes the square root of 
the mean square radius of gyration for the linear chaim 22 



FIG. 2: Double logarithmic plot of the form factor Pk{o) of 
an iV-noded Gaussian polygon under a topological condition 
K versus the variable u = qR,G,K- The curves colored with 
red, green, blue, magenta, and cyan correspond to the cases 
of all, 0, 3i, 4i and other, respectively. 



IV. SCATTERING FUNCTIONS 

For a random polygon with N nodes, we define the 
(single-chain) static structure factor g(q) by the Fourier 
transform of the pair correlation function as follows^ 

1 N 

g{q) = J dre^ r g(r) = - £ (exp(7g • (R m - R n ))) 

m,n— 1 

(2) 

We also call it the scattering function. The scattering 
function g(q) depends only on q = \q\, and we denote 
it by g(q). For a Gaussian polygon under a topological 
condition K , we denote the static structure factor or the 
scattering function by the symbol gi((q)- 

Let us introduce the form factor Pk(q) for a Gaussian 
polygon under a topological condition K as follows 



P K {q) 



9K{q) 

9k(0) 



(3) 



We have Pxiq) = 9K{q)/N from We have evaluated 
the form factor Pk{q) for the five topological conditions. 
Let us define variable u by u — qRo,K- The double log- 
arithmic plot of the form factor Pk {q) versus u is shown 
in Fig. |3 Here we recall that the form factor P a ii{q) was 
evaluated analytically in terms of the Dawson integral. 1 
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FIG. 3: Plot of (qR G ,K) 2 Pn{q) versus u = qR G ,K- The 
Kratky plots of the trivial, trefoil and figure-eight knots are 
represented by the green, blue and magenta curves, respec- 
tively. The red curve corresponds to the case of no topological 
constraint. 

In the region from u = Q up to tt = 2 or 3, the form 
factors Pk(q) for the five topological conditions overlap 
each other. For u > 5, the graphs of the different topo- 
logical conditions make parallel lines. The gradient is 
almost given by —2, which is consistent with the Gaus- 
sian asymptotic behavior. 

Let us discuss the Kratky plots of the form factors 
Pk{ci) for some different topological conditions. The 
plots of (qR~G,K) 2 PK{q) versus the variable u — qRc,K 
are shown in Fig. [2] Here we have numerically evalu- 
ated the Fourier transformations of the radial distribu- 
tion functions by an interpolation method^. 



K 


^pcak 


peak height 





2.08 


1.266 


3i 


2.09 


1.272 


4i 


2.09 


1.275 


other- 


2.09 


1.300 


all 


2.08 


1.287 



TABLE III: Peak position u pca k of the Kratky plot 
(<}Rg,k) 2 Pk{(£) versus u = qR G ,K for the Gaussian random 
polygon under a given topological condition K. The q value 
is given by an integral multiple of 0.01. 



plots are given by the same value of u for all the five 
topological conditions. The estimates of the peak posi- 
tions are given in Table ITTT1 The peak height depends on 
topological conditions. The Kratky plot for the trivial 
knot has the smallest peak height. The peak height for 
the trefoil knot is a little larger than that of the trivial 
knot. However, for the Kratky plots of the trefoil and 
figure-eight knots, the peak heights are given by almost 
the same value. 

The Kratky plots of Fig. [2] are not in contradiction 
with those of previous studies, and even generalize them. 
For lattice random polygons with N = 160, the Kratky 
plots were numerically evaluated for all polygons and 
knotted polygons, respectively^ Since the majority of 
polygons with nontrivial knots for N = 200 are given 
by those of the trefoil knot, the Kratky plots of all poly- 
gons and knotted polygons in Ref. approximately corre- 
spond to those of no topological constraint and the trefoil 
knot shown in Fig. [21 respectively. 

For u > 5, we observe that each of the Kratky plots 
of Fig. [21 approach constant values. We thus suggest an 
asymptotic behavior that for any topological constraint 
K the form factor Pk(m) should become close to that 
of the Gaussian polygon for u 3> 1, such as Pi((q) oc 
l/(qRG,aii) 2 - Here we recall that when u ^> 1 the form 
factor Pii n (g) of a Gaussian linear chain is approximated 
by Pun{q) « 2/{q 2 R G Hn ). We thus have 

R 2 

(qR GlK ) 2 P K {q) oc (4) 

U G,all 

The asymptotic constant value of the Kratky plot of a 
knot K should become smaller when the knot K becomes 
more complex. It is consistent with the observation of 
Fig. [21 that the order of the asymptotic constant values 
are the same with that of the values (R G K ) in Table [H 
for the three knots. For u > 10, systematic errors of the 
Kratky plots may be larger than statistical errors of the 
gyration radius Rg,k- 

The distribution function between two given nodes of 
a random polygon with a fixed knot is recently evaluated 
through simulation, and it is found to be close to the 
Gaussian onei^i The interpretation (@J should be also 
consistent with the observation. 



V. CONCLUSION 



For the graphs shown in Fig. [21 we observe that the 
mean square radius of gyration, (R G K ), plays a central 
role in the scattering function of the Gaussian polygon 
under a topological condition, K. In the small u region 
such as u < 2, the plots for the different topological con- 
ditions overlap completely, while for u > 2 the graphs 
become separate. Here, the peak positions of the Kratky 



We have evaluated the scattering functions of Gaussian 
random polygons with N — 200 under different topolog- 
ical conditions K. The Kratky plots have been obtained 
for the different K. They overlap up to u — 2, and they 
become separate for u > 2, and they approach constant 
values for u 3> 1. Several characteristic properties are 
explained in terms of the different values of the gyration 
radius Rg,k- 
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